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We give here a new third post-Newtonisn (3PN) spin-spin contribution (in the PN parameter 
e) to the accumulated orbital phase of a compact binary, arising from the spin-orbit precessional 
motion of the spins. In the equal mass case this contribution vanishes, but LISA sources of merging 
supermassive binary black holes have typically a mass ratio of 1:10. For such non-equal masses this 
3PN correction is periodic in time, with period approximately e _1 times larger than the period of 
gravitational waves. We derive a renormalized and simpler expression of the spin-spin coefficient at 
CN . 2PN, as an average over the time-scale of this period of the combined 2PN and 3PN contribution. 

5^ ' We also find that for LISA sources the quadrupole-monopole contribution to the phase dominates 

Oh' over the spin-spin contribution, while the self-spin contribution is negligible even for the dominant 

, spin. Finally we define a renormalized total spin coefficient a to be employed in the search for 

gravitational waves emitted by LISA sources. 
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I. INTRODUCTION 



Gravitational waves are generated by time- varying quadrupolar (or higher multipolar) mass (energy) configurations 
in general relativity. Although a convincing indirect proof for the existence of the theoretically predicted gravitational 
i ' waves is provided by binary pulsars [l[ , the most well-known being the Hulse- Taylor pulsar Q (others are also known 
0>[3|) a major effort is undergoing nowadays for the first direct detection of gravitational waves by the most sensitive 
Earth-based interferometric gravitational wave detectors LIGO @ and VIRGO @. 

These advanced instruments are able to detect in the frequency range 10—1000 Hz. Among their most likely sources 
are compact binary systems composed of neutron stars / stellar size black holes. In order to simplify the route towards 
the first detection, some of the relevant astrophysical details of the system are commonly omitted, like the spins and 
quadrupole moments. The detection effort concentrates on the signal emitted by point masses on quasi-circular orbit. 
With this, only the mass ratio should be monitored, while the inclusion of the astrophysical characteristics would 
induce many new parameters (additional six, for example, by the two spins). However the ambitious future aim of 
a gravitational wave based astrophysics would certainly require to include the spins and quadrupole moments of the 
| binary components in the data analysis. 

Supermassive black holes (SMBHs) of mass 3 x 10 6 -f- 3 x 10 9 solar masses are known to exist in the center of 
the majority of galaxies 0], p|, @, including ours. When galaxies merge, and the separation of the SMBHs is of 
a few parsecs, the SMBH binary enters a dissipative regime dominated by dynamical friction (an interaction of the 
SMBHs with the already merged stellar environment). As the relative separation further shrinks, this is followed by 
the gravitational radiation dominated dissipative regime, starting at approximately one hundredth of a parsec. 

During dynamical friction some of the orbital angular momentum of the binary black hole system is transferred 
to the stellar environment. In the process the stellar population at the poles of the system is ejected and a torus 
is formed [l(|, pd| . There had been a major worry, that the loss-cone mechanism for feeding stars into orbits that 
intersect the binary black holes is too slow and the SMBHs will stall in their approach to each other (the last parsec 
problem; for a review see fl2j). However recently more than one mechanism was proposed, which lead to a further 
approach of the SMBHs: (1) direct interaction with the surrounding stars slightly further outside [ll|; (2) relaxation 
processes due to cloud / star - star interactions [l3[, both repopulating the stellar orbits in the center of the galaxy. 
In a series of papers it has been recently shown that (3) even in the absence of two-body relaxation or gas dynamical 
processes, unequal mass and / or eccentric binaries with the mass larger than 10 5 solar masses can shrink to the 
gravitational wave emission regime of one hundredth of a parsec in less than a Hubble time due to the binary orbital 
decay by three-body interactions in the gravitationally-bound stellar cusps [14j . Alternatively, by considering (4) the 
three accretion disks: one around each black hole and a third one, which is circumbinary, it has been shown in (l5j 
that the balance of the orbital angular momentum changing processes (the circumbinary disk removes orbital angular 
momentum from the binary via the binary-disk resonant interaction, however the mass transfer to each individual 
black hole adds orbital angular momentum to the binary) is such that these binaries will merge within a Hubble time 
driven by this mechanism. 
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In consequence by one mechanism or another the SMBHs would be able to approach each other to distances smaller 
than approximately one hundredth parsecs, when the gravitational radiation becomes the dominant dissipative effect. 
The gravitational waves in the frequency range 10~ 4 — 10 _1 Hz emitted by such a system of coalescing SMBHs will 
be captured by the forthcoming LISA space mission (l6j . 

There are a few important comments to make on these SMBH LISA sources. 

(al) The mass distribution of the galactic central SMBHs shows a broken powerlaw [13], [Hjl, confirmed by most 
recent work [l^ | and a recent observational survey [20l | . 

(a2) Let v — m^jmi be the mass ratio of the binary with masses m-i < mi. Thus < v < 1. Next we define 
m — mi + rri2 as the total mass, and = m\milra as the reduced mass, ranging between v for v — > to 1/4 for equal 
masses. Further we introduce the symmetric mass ratio r\ — fx/m = (2 + v + J^ -1 ) (for small mass ratios w — > v). 
The probability for a specific mass ratio for SMBH encounters can be derived as an integral over the black hole mass 
distribution, folded with the rate to actually merge (proportional to the capture cross section and the relative velocity 
for two galaxies). By this method it has been estimated in [2l[ that the most likely mass ratios are in the range 
v 6 (1/30, 1/3). Therefore the most typical mass ratio for LISA sources would be v — 10 _1 , to be considered in this 
paper. 

(bl) The spin of individual black holes has a quite high value. As the surrounding matter is accreted into the hole, 
the angular moment transfer spins up the black hole, close to the maximally allowed value. By taking into account 
the angular momentum transfer from accreting matter alone, the dimensionless spin parameter \ = S/ra 2 (in units 
c = 1= G) grows to the maximally allowed value of 1 even if the initial state of the black hole was non- rotating [2l| • 
By taking into account the energy input of the in-falling (horizon-crossing) photons emitted from the steady-state 
thin accretion disk, a torque counteracting to the one due to mass accretion has been found [23I ] , reducing the limiting 
value of the dimensionless spin parameter to 0.9982. The description of the process can be further enhanced by (i) 
including the contribution of the open or closed magnetic field lines in the vicinity of the disk [241, (ii) discussing 
a symbiotic system (introduced in (25j) of a rotating black hole, a magnetized thin accretion disk and jets in the 
magnetosphere of the hole (26J. With any of these refinements, the end state of the black hole will stay very close to 
the extremal Kerr limit. 

(b2) In consequence the spin of the individual components of a binary black hole system, which is formed typically 
by capture events (even for stellar mass black holes the emergence of binary systems by stellar evolution is quite 
unlikely) will be also high. 1 

We conclude from these that any reasonable analysis of the gravitational waves produced by central galactic SMBH 
mergers should take into account both the individual spins Si of the black holes and their typical mass ratio. 

In the inspiral phase of the merger, which can be described by post-Newtonian (PN) techniques and lasts until the 
innermost stable orbit (ISO) is reached, the spins and quadrupole moments induce a quasi-precessional evolution of 
the spin vectors [3lj |. The leading contributions come from the SO precessions at 1.5PN order. The interaction of the 
individual spins lead to a 2PN spin-spin (SS) quasi-precessional contribution to the dynamics. Similar 2PN quasi- 
precessional contributions arise, when one binary component is regarded as a monopole moving in the quadrupolar 
field of the other (quadrupole- monopole, QM contribution). This angular evolution of the spins is decoupled from 
the radial motion, which was completely integrated in terms of a generalized Kepler equation [32| . The inclusion 
of the angular evolution (spin precessions) into the description besides increasing the accuracy, also breaks many 
degeneracies among parameters, greatly improving (relative to previous estimates) how well LISA will be able to 



measure the masses, the locations of binaries on the sky, and the distance to the source |33| . The gravitational 
waveforms for spinning compact binaries with precession effects through 1.5PN order were also recently analyzed 
using spin- weighted -2 spherical harmonics, and applied to spinning, non-precessing binaries [34| . 

A careful and detailed discussion, following the assumptions of Ref. (35[, by taking into account both the SO 
precession and the dissipative effect of gravitational radiation (averaged over one radial period) for the above typical 
mass ratio has already resulted in two achievements [2l| . The first is a quantitative description of the spin-flip process, 
visible in the X-shaped radio galaxies. Secondly, the precessional phase of the merger of two black holes, occurring 
prior to the spin-flip, was suggested to be the source of the superdisk, visible in radio galaxies. The precessing jet 
acts as a super-wind separating the two radio lobes in the final stages of the merger. According to this model such 
radio galaxies are candidates for subsequent SMBH mergers [2l[ . 



How large is the final spin after the merger occurs, is still the object of debate. Recent numerical work on the merger of two black 
holes with spin has predicted (beside recoil effects arising for particular configurations [27|, interesting in their own), the final spin state 
of the merged black hole. Extrapolating from the particular configurations analyzed, an analytical formula was also proposed for the 
final spin l28ll . This formula has been employed in the analysis of the final spin as arising from a combination of accretion and merger 
sequence |29( . Alternative analyses of the final spin in mergers have been performed by analytical or numerical methods, yielding various 
predictions for the values of the final spin [30l | . 
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The actual detection of gravitational waves requires the precise knowledge of the waveforms. This is determined 
both by the relative geometry of the sources and detectors, and by the evolution of the sources. In the latter, the 
accumulated orbital phase <f> plays a central role, as any effect that causes the template to differ from the actual 
signal by one cycle over the 500 to 16,000 accumulated cycles in the sensitive bandwidth will result in a substantial 
reduction in the signal-to-noise ratio [361 ] . Therefore the phase of the gravitational waves emitted by a compact binary 
is required to high PN accuracy. Various contributions up to 3.5 PN orders are known. Among these, there are 
important spin and mass quadrupole contributions, like SO at 1.5PN, also SS, QM and SS-self contributions at 2PN. 
(The latter arise from the self- interaction of the spins [37j .) We note here, that alternative definitions of the phase 
in the spinning case were introduced (the orbital phase with respect to the ascending node [38|) and the rigorous 
waveform containing SO contributions has been also computed [391 ]. 

Besides possible quadrupolar irregularities in the mass distribution of the binary components (which are however 
quite unlikely due to the strong self-gravity of these objects yieldin g ra pidly to a spherical configuration), and the 
mutual quadrupolar deformation of them (shown to be quite small in [40j] ) , the proper rotation is the effect generating 
a considerable quadrupole moment for each black hole. The quadrupolar deformation can be characterized by the 
coefficients pi = Qi/rriim 2 , with Qi the quadrupole scalar. If the quadrupole moment originates entirely in the rotation 
of the compact object (what we shall assume in what follows), then Qi ~ — axfraf , with the parameter a £ (4, 8) for 
neutron stars, depending on their equation of state. Stiffer equations of state give larger values of a, while a = 1 for 
rotating black holes [4l[. Here Xi = Si/ m i are the dimensionless rotation parameters and the negative sign is because 
the rotating compact object is centrifugally flattened, becoming an oblate spheroid. 

In what follows, we will consider maximally rotating SMBHs (xi ~ 1; thus Si mf), with their quadrupole moment 
entirely determined by the rotation. Therefore pi — —Sf/mfm 2 s» — (rrii/m) , thus 

pi » -{1 + vT 2 , (1) 

P2 » -(1 + ^-T 2 - (2) 

The ratio of the spins of the compact objects (characterized by Ri « mi and evaluated for comparable and high 
rotation velocities Vi ~ V% « 1 according to our assumption of maximal rotation) can be expressed as 

m 2 \ _ 2 



, , , ^, (3) 
while in the same approximation the ratio of the spins with the magnitude L of orbital angular momentum L becomes 

^ - ^=nfiv 2 -^ v v (4) 



L fxrv \ r ) \v J mi 

e 1 ^- 1 , (5) 



Si 5*2 <S*i 

T ~ ~ls~ 2 



where we have employed the definition of the PN parameter e rj m/r w v 2 . As remarked in [2l| . in case of a small 
v « r\ the ratio S-z/L is shifted towards even smaller values (therefore S2 <C L during all stages of the inspiral), 
however the ratio of the spin of the dominant compact object to the magnitude of the orbital angular momentum can 
change drastically. Indeed, it is determined by the relative magnitude of the small parameters e and v. As e increases 
during the inspiral, whenever v falls in the range of e 1 / 2 , an S\ « L epoch is reached, which follows the initial epoch 
with Si < L and precedes the forthcoming Si > L epoch. This happens for the mass ratios v £ (1/30, 1/3). 

We also note that in previous works we have found useful to choose as some of the parameters for the new spin 
degrees of freedom the relative angles Ki between the spins and orbital angular momentum and the angle 7 between 
the spins. This is because both the quasi-precessional equations for the spins and the equations describing the 
gravitational radiation losses up to 4.5 PN orders close in these variables, the energy E and magnitude of the orbital 
angular momentum (averaged over one radial orbit) L [42j , [43| , [3 , [45| . 

In this paper we will review the contributions induced by the spins and quadrupole moments to the accumulated 
orbital phase of the binary and we propose modifications in this description. In order to do so, in Section II we 
review the existing contributions up to 3.5PN orders. Then we represent graphically the dependence of the QM 
contribution on the parameter space (ki, k^). We illustrate that due to the involved numerical coefficients, the self- 
spin contribution, although of 2PN order, is much smaller, than the QM contribution, thus it can be safely dropped 
at this accuracy. Both the QM and the SS-self contributions are discussed for either equal or unequal masses. The 
dependence of the SS contribution on the parameter space (ki, K2) and on the relative spin configuration, expressed 
in terms of the relative spin azimuthal angle Aip is discussed only in the equal mass case here. (Due to a spherical 
triangle identity this variable is connected to the previously mentioned angles Ki and 7.) 
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In Section III we discuss the consequences of the leading order precession effects on these coefficients. We show that 
due to the SO precession the SS contribution to the phase receives 3PN time-dependent corrections, which happen 
to cancel in the equal mass case. At non-equal masses however these contributions vary periodically, with a period 
T3PNSS which is e~ x times longer than the period of the gravitational wave T wave . 

Then in Section IV we look for LISA sources of typical mass ratio v = 1CP 1 . We represent here the detailed time 
evolution of the total spin-spin coefficient on the time scale Tspnss in the parameter space (/ei, K2). Next we repeat 
this for the sum of SS and QM contributions, showing that the QM effects dominate and the SS effects modulate. 
We are then able to define a renormalized SS coefficient at 2PN, which can be considered constant up to 3PN on a 
time-scale T^p^ss = £~ l T waV e- This coefficient is simpler than its unnormalized counterpart. 

In Section V we compare the number of cycles left for typical LISA sources, considering the following cases: (a) 
the 2PN accurate spin-spin coefficient, previously employed in the literature, with the relative spin angles 7, K{ being 
constants at this accuracy (the precessions contributing only with higher order terms, thus dropped here); (b) our 
renormalized spin-spin coefficient (constant up to 3PN), without the precession equations, which arc already taken 
care of to this accuracy by our renormalization procedure; (c) the numerical evolution, with time-varying spin-spin 
coefficient, found by including the precessions. 

Finally in the Concluding Remarks we summarize our findings and define the renormalized value of the total 
coefficient encompassing both the QM and SS contributions, which due to its simplicity, time- independence up to 
3PN, and much better agreement with the numerical results is more suitable to be used in future LISA data analysis 
than the presently available expression. 

Throughout the paper we use units G = 1 = c. 



II. THE PHASE OF GRAVITATIONAL WAVES TO 3.5 PN ACCURACY 

The integrated orbital phase of gravitational waves has been computed to high accuracy, with contributions both 
from general relativistic corrections to the Keplerian motion and other Newtonian contributions due to the finite-size 
of the compact binary components. The respective contributions were recently summarized in [4r| . [37| . and (47j . Up 
to the desired 3.5 PN accuracy the phase can be formally decomposed as: 



<t> = <t>c + 4>N + <t>lPN + 4>l.hPN + <t>2PN + <p2.5PN + 4>3PN + 03.5PJV 

where <j) c is a constant, and the lower order contributions are completely known as: 

1 
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(7) 
(8) 
(9) 
(10) 
(11) 



Here the dimensionless time parameter r = r/(t c — t)/5m is decreasing until the merger. The contribution at 1PN is 
purely relativistic, computed in 14811 . At 1.5 PN appear the leading-order tail contributions [49j-[52j, together with 
the spin-orbit contributions [54j |. [55j |. encompassed in the coefficient 
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At 2PN there are both relativistic contributions [53| and finite size effects, collected into the coefficient 



a SiS 2 



CTSS-self 



a SiS 2 + a SS-self + &QM , 
S S 

t(— 247cos7 + 721 coski COSK2) 

48r/m 4 



static- 
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VQM = ~ ~ ^2 Pi (3 COS 2 Ki - l) 



(13) 
(14) 

(15) 
(16) 



551 ] . self-spin [37| and quadrupole-monopole [4l| interactions. 
if on the relative angles Ki is represented both for equal 



The enlisted contributions originate in the spin-spin [54 

The dependence of the spin coefficients <jqm and ass- 
masses and v = 10 _1 on Fig. [TJ under the assumptions of maximal rotation and quadrupole moment generated purely 
by rotation. Although formally both contributions are of 2PN order, the numerical values of the coefficients involved 
imply that the quadrupole-monopole contribution clearly dominates over the self-spin contribution, which only takes 
values in the narrow interval <Jss-seif € (1/32, 7/96) 52i=i 2 ( m i/ m ) ■ 

The coefficient <tqm strongly depends on the spin configuration, in particular in the v = 10 _1 case only on the 
relative angle between the dominant spin and orbital angular momentum. By comparing the two panels of the figure, 
we see that the <jqm contribution increases considerably with decreasing mass ratio. 



V = 1 



10-1 





FIG. 1: (Color online). The quadrupole-monopole contribution ctqm (blue wavy surface) represented for equal masses and 
v — 10" 1 in the parameter space k%). The coefficient oqm strongly depends on the spin configuration, in particular in the 
v — 1CP 1 case only on the relative angle between the dominant spin and orbital angular momentum. The oqm contribution 
increases considerably with decreasing mass ratio. For comparison, the much smaller self-spin contribution ass-seif is also 
represented (green levelled surface). The values of ass-seif are positive, but well below 1. The plots are for maximally spinning 
binaries, with the quadrupole moment induced by rotation. 



In the equal mass case we also show the dependence of ts 1 s 2 an( i °f the total a on the parameter space («i, K2) 
for various relative spin azimuthal angles (Figs. [2] and [3]). Here = V>2 — 4>i is the difference between the 
azimuthal angles ipi of the spin vectors. As the vertical scale for all these figures is the same, we can clearly see 
that the quadrupole-monopole contribution dominates for a huge part of the parameter space. Also, in general a ^ 
and its positive values are higher than the negative ones due to the QM contribution. However there is a small set of 
particular configurations for which the contributions to a cancel out each other. 

The coefficients crs 1 s 2 an d cr in the non-equal mass case will be discussed in the remaining Sections of the paper. 
They are more complicated because (anticipating the results of the next Section) the relative spin azimuthal angle 
Atp becomes time-dependent for v ^ 1 due to precessions. 

Finally at 2.5 PN there are contributions to <fi from the first PN correction to SO [56j], [47|, expressed by 

2 r, 
V—v Oi COS Ki 
PPN = > , 2 
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FIG. 2: (Color online). The dependence of the spin-spin parameter (JSi S2 ( re d wavy surfaces) on the angles Ki for equal masses 
and various values of the relativ spin azimuthal angle Aip = 0, 7r/2, tt, Sn/2. For comparison the zero level surface is also 
represented. Plot for the maximal rotation case. 




FIG. 3: (Color online). Same as on Fig [2] but for the total coefficient a (lilac wavy surfaces). Although there are particular 
configurations (at the intersection with the zero level surfaces) for which the various contributions to a cancel out each other, 
for the vaste majority of the parameter space a 7^ 0. The positive values are typically higher than the negative ones due to the 
QM contribution. 



and tail effects [571 ] . We note that the logarithmic terms arising from the latter can be eliminated through a modified 
definition of the phase 58] 

The list of contributions known at 3PN and 3.5PN are not exhaustive. At 3PN there are known relativistic and 
tail contributions HfJ 
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where C = linin^oo Efc=i 0-/^) — l nn ] ~ 0.577 is the Euler-Mascheroni constant 2 . However the contributions from 
the PN corrections of the SS, self and QM contributions are not known. 
At 3.5PN only the tail contribution has been calculated [46| 



1 / 188516689 488825 141769 



'J3.5PN — 



r) - ^-^if ^" 1/4 , (19) 



j] V 173408256 516096' 516096 

but clearly there are numerous other contributions, as the 2PN correction to the SO contribution, SO-SS, SO-QM 
couplings etc. 

It is the main purpose of this paper to add new spin-spin contributions at 3PN, arising from the spin-orbit precession 
of the spin vectors and discuss their implication to the definition of an effective a at 2PN. 



III. LEADING ORDER PRECESSIONAL CONTRIBUTION TO a 

The s pin vectors undergo a complicated quasi-precessional motion [3l| , and the angles Ki and 7 evolve accordingly. 
In Ref. [42j it was derived that the leading order contributions to the quasi-precessional variation of the angles Ki is 
of e 3 / 2 order, while 7 evolves in the e order. Supplementing the SO and SS contributions already considered in Ref. 
[42I ] with the QM contributions (with the quadrupole moment arising from pure rotation), by employing Eqs. ([I])-© 
and the method described in a footnote of Ref. [44j |. according to which O(Sx) — 0(i)re -1 / 2 , we have checked that 
these estimates are not changed, although the leading order changes in Ki acquire QM type contributions. However 
the leading order change in 7 arises solely from the SO precession. 

With the angles k, and 7 non-constant, the coefficients j3, a and (3pn also evolve in time. From among the relative 
spin angles the contributions uqm and <rss~seif contain only Ki. As these start to vary at e 3 / 2 order, the variations 
of Ki would contribute only at 3.5PN orders. The same is true for the Ki-dependence of osiS 2 - However the variation 
at e order of 7 induces 3PN corrections to ersiS 2 - Thus the leading order precessional contribution to a is caused by 
the leading order precessional evolution of the spin angle 7, of SO origin. 

Therefore in what follows we derive the time dependence of 7. From the SO precession of the spin vectors (the SO 
contribution in Eq. (2.17c) of Ref. [42j], with the change of notation 77 — > v) we find 



(cos 7)'= o sin k\ sin K2 sin Aip . (20) 



3L {y 
~2r 

By taking the time derivative of the spherical cosine identity 

cos 7 = cos Ki cos Ki + cos Aip sin Ki sin k 2 ■ (21 

we get 



sm7 — = sin K\ sin K2 sin Aip— Alp 
at ~ dt 



- cos Atp 



d d 

COS «i — (COS K2 ) + COS K2 — (COS K\ } 

(Jib (JjV 



(22) 



The second term can be dropped as we need to consider this relation only to e accuracy. A comparison of Eqs. ([20 
and (22]) gives: 



whenever Kj, ^ 0. For equal masses Aip =const, therefore 7 is also unchanged by precessional effects to e order. 
Therefore <Js 1 S 2 can b e considered constant even at 3PN for v = 1. 

For any v ^ 1 Eq. (|23[) can be integrated in the following way. First we note that in order to calculate the leading 
order change in Atp, all quantities are needed only to Keplerian order. We also specialize to circular orbits r = a, as 



2 The Euler-Mascheroni constant appears, among other places, in expressions involving the exponential integral, the Laplace transform 
of the natural logarithm, an inequality for Euler's totient function and solution of the second kind to Bessel's equation [59ll . 
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required by the derived expression of the phase. We then pass to integration in terms of the true anomaly x, defined 
as 

^ = — (24) 
dt ~ fxa 2 ' 1 > 

and employ the Kepler equation to establish the connection between the time and the eccentric anomaly £ (which to 
leading order and circular orbits coincides with the true anomaly x): 

n(t-t )= X , (25) 
where n — 2ir/T or bit- By choosing to = 0, we obtain the required explicit time dependence of 7: 

3/in (y* 1 — v) 



COS 7 = COS Ki COS K2 + COS 



(AV)r 



2a 



sin«i sinK2 , (26) 



which gives a new time-dependent 3PN spin contribution to <7SiS 2 - 

The angle 7 (and the respective 3PN contribution to <7s 1 s 2 ) nas a periodicity with period 

iwa 2a 

T3PNSS = 1 _, r~ = , _, zT orb it ■ (27) 

3/i [y 1 — v) n 3/i (y 1 — v) 

With the approximation valid for quasicircular orbits: T or ut — 2T TOOt)e , also employing e « raja and t^ -1 = m//i = 
(2 + ^ + ^ _1 ), we get: 



r 3 Pivss 4 4 2 + v + v- 



T waV e 3 (v 1 - V) ET] 3 (l/ 1 - V) S 



(28) 



It has been shown in [21(, that the typical mass ratio at galactic black hole encounters is between 1/3 to 1/30, such 
that the mass ratio v = 10 _1 can be considered fairly representative. We adopt this value in what follows. Then 

^™ * 1.63s- 1 . (29) 

We have therefore found that for the representative mass ratio the time-scale of variation of the angle 7, (and the 
related change in <Js 1 s 2 ) is °f £ times larger order, than the period of gravitational waves. 



IV. THE TIME-EVOLVING a Sl s 2 AND a IN THE UNEQUAL MASS CASE v = 10 _1 
We remark on Eq. (|26[) that on a time-scale Tzpnss-, the time-dependent contribution averages out: 



COS 7 = COS K\ COS K2 , 

therefore the average of the spin coefficient (Js 1 s 2 becomes 

a S!S 2 = -3 4- COS Ki COS K 2 



(30) 



(31) 



We call this the renormalized spin-spin coefficient in the phase of gravitational waves. A comparison with the standard 
expression (|14p shows that the renormalized expression is simpler than the original one. 

The spin coefficient (?SiS 2 to 3PN accuracy can be then rewritten as a sum of the renormalized coefficient and a 
time- varying part: 



<T S 1 S 2 = CTSiS 2 + 5(T , 

247S'iS*2 



5a = 



487/TO 4 



■ sin K\ sin K2 cos 



, . . 3un (v 1 — v) 
2a - * 



(32) 
(33) 



The second term is periodic (with period T 3 p^gg w e l T wave for the typical binaries with mass ratio v — 10 1 ), and 
as such, it averages out over a time scale e~ l T wave . 
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v = 1CH 




FIG. 4: (Color online). The time-evolving spin-spin parameter as 1 s 2 for mass ratio v = 10 1 represented every quarter period 
during T 3P xss- 



v = 10-1 




FIG. 5: (Color online). The time-evolving parameter a for mass ratio v — 10 1 represented every quarter period during 
Tspnss • One can see that the dominant contribution is from the quadrupole-monopole interaction. 

On Figs S] and [5] we have represented the time-evolving parameters and o for the mass ratio v = 10 _1 in a 

time sequence sampled every quarter period during T^pnss- 3 While the importance of the spin-spin contribution 
decreased as compared to the equal mass case (the scales on Figs [2] and 0] are the same for easy comparison), the 
quadrupole-monopole contribution is enhanced by a small v (Fig[T]). As a result, for LISA sources the dominant con- 
tribution comes from the quadrupole-monopole interaction and the time-varying spin-spin contribution only induces 
a wobbling of the crests determined by the quadrupole-monopole interaction (compare Figs [1] and [5]) . 



3 Animated gif-s showing the temporal evolution of both cs 1 g 2 and a are available at [60l |. 
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V. ANALYTICAL AND NUMERICAL CONTRIBUTIONS TO THE ACCUMULATED NUMBER OF 

GRAVITATIONAL WAVE CYCLES 



The phase of the gravitational waves ((6|) emerges by two integrations over time from the expression of the time 
derivative of the orbital angular frequency Co, see Ref. [37j ■ The spin coefficient as t s 2 nf st appears in the latter 
expression and to a 2PN accuracy it can be considered constant, hence it also emerges in the phase. 

To a higher accuracy the coefficient acquires time-dependence from the precession equations, which should change 
the results of the integration. Alternatively, the contribution of the leading order SO precession can be taken into 
account by the averaging procedure described in the previous section. Then the spin coefficient as[ s 2 , regarded as 
a constant, formally leads to the same expression of the phase and accumulated number of gravitational wave cycles 
Af = (4> c — 4>) /it as Eq. (12) of Ref. [33], with the renormalized coefficient asiS 2 replacing asiS 2 - Obviously this 
result is approximate rather then exact. The question comes, how good is the approximation? 

As the exact time integrations turn out to be cumbersome and the result is not expressible in terms of elementary 
functions, in order to answer this question we calculate numerically the corresponding contribution to the accumulated 
number of gravitational wave cycles. We proceed in the following way. Eq. (7) of Ref. [33] gives the evolution of 
the orbital angular frequency to in time to 2PN accuracy. By inserting the time dependent expression (|3"2")l into this 
equation, and keeping only the Newtonian and the time-dependent correction of the spin-spin terms, we obtain: 



d(muj) ii /3 5 

— = —96 mw ' — 96 (moo) da 

dr 



The solution of this equation can be also split into the known Newtonian contribution and a term of 5a order: 

r -3/8 



■ mivsa (t) 



(34) 



(35) 



Here the dimensionless time t is defined as t = rj (t c — t) /5m, with (t — t c ) the time until the coalescence. (Note that 
we are discussing only the precessional corrections to the spin-spin contribution in this paper, without attempting to 
encompass other cumbersome post-Newtonian corrections of the spin-spin effects arising from the dynamics.) 

Rather than giving the number of cycles left in terms of r, we express it as an integral over the orbital angular 
frequency 



N = M N + Nsa = - u (t) dr = - / —doj . 

7T .L. 7T ./,.,. LO 



(36) 



Here Co is given by Eq. (|34j) . In order to express 5a (to), from Eq. (|33|) by employing Kepler's third law to 2 a 3, = m 
and n = lo (valid to leading order), we obtain 



4nm 



3PNSS 



3?7 (y 1 — v) (mLo) 5 ^ 
We also express t as function of lo to leading order from Eq. ([35]) as 



t-t c 



5m 



t \-8/3 



Thus 



247SiS 2 . 

oa — — sin Ki sin K2 cos 

48ryTO 4 



/a I \ , 3? 7 (y 1 - v) tc , ,5/3 

(AV>) H —^z — — {mu)) 



2m 



29 ( TOW ) 



(37) 



(38) 



(39) 



Inserting the expression (|34|1 . with 5a given by Eq. (f39|) into Eq. (|36|) we obtain the number of cycles in the frequency 
range (oji,ujf). This enables us to give accurate numerical results for the accumulated number of cycles. 

For comparison we select three LISA sources with unequal mass given in Tables I and II of Ref. [62| and we enlist 
in Table U the corresponding spin-spin contributions calculated for particular (constant) values of the angles Ki in 
three different ways: (a) by employing Eq. (10) of Ref. (37j : (b) in the same way, but with our renormalized spin-spin 
coefficient as 1 s 2 replacing er^Sj,; (c) employing a numerical integration based on the exact method described in this 
section. For comparison we also give the corresponding Newtonian, PN and QM contributions. The data confirms 
that for these non-equal mass binaries the QM contribution dominates over the spin-spin contribution. It is also 
immediate to see that the renormalized spin-spin coefficient leads to much closer values to the results of the numerical 
integration, then the 2PN-accurate coefficient. 
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TABLE I: Newtonian, PN, QM and spin-spin contributions to the accumulated number of gravitational wave cycles for three 
LISA sources with unequal mass, for particular values of the angles Ki and (A^) Q . The relevant frequency ranges (with 
/ = ijj/2n) are also indicated. The spin-spin contributions are computed in three distinct ways (a) in a 2PN rigurous formalism 
(A/"o- ss ); (b) with the renormalized spin-spin coefficient (A/"a SS ); and (c) taking into account the exact precessional evolutions 
up to 3PN by numerical integration (Ma- ss +sa- S s( t ))- The numbers obtained with the renormalized coefficient are within 0.02 
cycles agreement with the numerical data, presenting an improvement up to 0.6 cycles over the 2PN accurate expression. The 
numbers in the table confirm that for unequal masses the QM contribution dominates over the spin-spin contribution. 



(AV) = 15° 


1.4M© 


10 4 M Q 


10 6 M© 


Kl = K2 = 30° 


400M Q 


10 5 M Q 


1O 7 M 


MHz) 


4.601 x 10 -2 


4.199 x 10 -4 


2.361 x 10 -5 


ffin(Hz) 


1.000 


3.997 x 10" 2 


3.997 x 10" 4 


Mn 


2294624.602 


21055.705 


1174.257 


Mpn 


35363.312 


677.247 


114.536 


Ao-Qjvf 


-432.606 


-13.712 


-6.028 




-2.986 


-2.678 


-1.177 


Mg ss 


-3.588 


-3.218 


-1.415 




-3.588 


-3.207 


-1.404 



VI. CONCLUDING REMARKS 



We have discussed the contributions induced by the spin-spin and quadrupolc-monopole coupling to the accumulated 
orbital phase. We have shown that the contributions from the self-coupling of the spins are much smaller than the QM 
and SS contributions in all cases (FigQ}, therefore they can be dropped at a 2PN accuracy. The QM contributions 
depend on the parameter space (ki, K2), while the SS contributions on (ki, K2, At/;) [or alternatively on (ki, K2, 7)], 
and were represented on these parameter spaces. 

For equal mass compact binaries we have shown that the QM contributions (left panel of Fig [1]) and the SS 
contributions (Fig[2]) are of comparable size, although the SS contributions are slightly higher. These two contributions 
(at least for positive values) happen to roughly amplify each other, thus their combined contribution (Fig [3]) has 
approximately twice as large positive amplitude, than the individual contributions. There is also a set of measure 
zero in the parameter space, for which the total contribution vanishes. 

For typical LISA sources with unequal masses v — 1CP 1 the SS contribution becomes time-dependent due to the 
variation of the angle 7 by SO precession. The time-evolution of the SS coefficient was shown on Fig0J By comparing 
with the QM contribution for v = 10 _1 (right panel of Fig [I]), we see that the SS contribution is much smaller. As 
consequence, when we add these contributions, the combined coefficient shows the pattern of the QM contribution, 
modulated in time by the SS contribution (Fig [SJ) . 

The time variation of the SS component was shown to be periodic, with a period e _1 times longer than the period 
of the gravitational wave. By separating from the SS contribution a time-varying component which averages to zero 
on this time-scale, we may call the rest as a renormalized SS coefficient as 1 s 2 - Besides of being constant up to 3PN, 
it has a simpler expression, than the previously known 2PN SS coefficient 05^ s 2 . 

Therefore we find useful to introduce in the unequal mass case the total renormalized coefficient W = uqm + <?s7s7+ 
ass-self for maximally rotating black holes with quadrupole moment due to rotation. It has the expression 



5 



3 cos 2 ki — 1 



„ 79 

2 H cos Ki cos Ko 

20 



6 COS 2 Ki 



(40) 



which was obtained by keeping only the leading and first order QM and SS contributions in v. The renormalized a is 
constant in time up to 3PN and as such is more convenient to use in wave detection, than a, which is constant only 
to 2PN. For all numerical examples studied, the renormalized coefficient gave results within 0.02 cycle accuracy as 
compared to the numerical results. This represented an improvement of 0.4-^0.6 cycles over the previously employed 
2PN accurate coefficient. 

The dependence of a on the parameter space («x, ^2) is shown on Fig [5] It is very similar to the pure quadrupole- 
monopole effect, represented in the second panel of Fig [I] which is slightly modulated by the renormalized spin-spin 
contribution. 
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i' = 10 _1 




FIG. 6: (Color online). The renormalized parameter a for mass ratio v = 10 _1 . The dominant contribution comes from the 
quadrupole-monopole interaction (compare with the second panel of Fig [TJ , which is slightly modulated by the renormalized 
spin-spin contribution. 

Due to the arguments presented in this paper we propose to replace a by a in the expression of the gravitational 
wave phase, Eq. (|10p , whenever it is applied for LISA sources. 

In a similar way to the time dependence and renormalization of a, a time dependence and renormalization of the 
spin-orbit contribution (3 can be derived, the calculation of which involves to integrate the combined SO, SS and QM 
quasi-precessional time evolution of the angles Ki at e 3 / 2 order. This has been done for the equal mass case [6l| . 
however it seems to be far from obvious for v ^ 1. Therefore we defer the study of this effect to a future work. 
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